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ABSTRACT:  Boron  carbide  (B4C)  is  a  hard  material  whose  value  for  extended 
engineering  applications  such  as  body  armor;  is  limited  by  its  brittleness  under  impact. 
To  improve  the  ductility  while  retaining  hardness,  we  used  density  functional  theory  to 
examine  modifying  B4C  ductility  through  microalloying.  We  found  that  replacing  the 
CBC  chain  in  B4C  with  Si— Si,  denoted  as  (BuCp)  — Si2,  dramatically  improves  the 
ductility,  allowing  a  continuous  shear  to  a  large  strain  of  0.802  (about  twice  of  B4C  failure 
strain)  without  brittle  failure.  Moreover,  (BnC)— Si2  retains  low  density  and  high 
hardness.  This  ductility  improvement  arises  because  the  Si— Si  linkages  enable  the 
icosahedra  accommodate  additional  shear  by  rotating  instead  of  breaking  bonds. 


Sio^° 


Substitute  C-B-C  chain  with  Si2  pair;  leading  to 
lower  density,  and  increased  ductility. 


SECTION:  Plasmonics,  Optical  Materials,  and  Hard  Matter 


Superhard  materials,  such  as  diamond,  cubic  boron  nitride, 
and  boron  carbide  (B4C),  exhibit  many  valuable  properties: 
high  melting  temperature,  large  compression  strength,  chemical 
inertness,  and  high  thermal  conductivity,  making  them  of 
practical  importance  for  science  and  engineering  applica¬ 
tions.1-6  However,  they  are  brittle  with  low  fracture  toughness. 
This  serious  flaw  prevents  many  engineering  applications,  in 
particular  involving  extreme  environments  of  high  stress  and 
hypervelocity  impact.  Efforts  to  improve  the  fracture  toughness 
include  designing  functional  microstructures  and  forming 
porous  structures.7  10  However,  improvement  of  fracture 
toughness  is  often  at  the  expense  of  hardness  and  wear 
resistance.  To  achieve  new  materials  for  extended  engineering 
applications,  we  want  to  redesign  hard  materials  to  be  ductile 
without  sacrificing  these  properties. 

Boron  carbide  (B4C)  is  extremely  hard  (Vickers  hardness  30 
GPa),  although  less  hard  than  diamond  (Vickers  hardness  115 
GPa)  and  cubic  boron  nitride  (Vickers  hardness  76  GPa). 
Because  of  its  novel  properties  of  high  hardness,  low  density, 
and  abrasion  resistance,  B4C  is  widely  used  in  refractory 
applications,  abrasive  powders,  and  body  armor  and  as  a 
neutron  radiation  absorbent.11  21  However,  B4C  fractures 
easily  just  above  the  Hugoniot  elastic  limit  (HEL,  ~20  GPa) 
at  hypervelocity  impact  and  high  pressures.22  25  This  low 
fracture  toughness  results  from  formation  of  local  amorphiza- 
tion  shear  bands,16'26-28  for  which  the  atomistic  mechanism  has 
been  elucidated  from  combining  experimental  and  computa¬ 
tional  studies.16  18,-7  30  In  particular,  we  examined  the 
response  of  B4C  to  shear  along  11  slip  systems  and  found 
that  the  slip  system  with  the  lowest  shear  strength,  (0111)/ 
<  1 101  >,  leads  to  a  unique  plastic  deformation  before  failure  in 
which  a  boron— carbon  bond  between  neighboring  icosahedral 
clusters  breaks  to  form  a  carbon  lone  pair  (Lewis  base)  on  the 
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C  within  the  icosahedron.30  Further  shear  then  leads  this  Lewis 
base  C  to  form  a  new  bond  with  the  Lewis  acidic  B  in  the 
middle  of  a  CBC  chain.  This  then  initiates  destruction  of  the 
icosahedron  and  subsequent  failure. 

We  take  the  hard  material  B4C  as  the  prototype  and  examine 
how  to  improve  its  ductility  through  micro  alloying.  Using  the 
PBE  flavor  of  density  functional  theory  (DFT),  we  examined 
the  changes  of  mechanical  properties  including  ductility, 
density,  and  hardness  by  modifying  the  CBC  chains  connecting 
the  icosahedra  in  B4C.  We  first  assessed  the  ductility  of  various 
modified  structures  using  the  empirical  criterion  that  the  ratio 
of  the  bulk  modulus  to  the  shear  modulus  (B/G)  should  be 
better  than  B4C  (>1.20).  Then,  we  selected  the  best  cases  with 
improved  ductility,  low  density,  and  high  hardness  and  sheared 
them  along  the  amorphous  slip  system  in  B4C  to  examine  the 
shear  failure  model  and  ductility  of  new  systems. 

The  description  of  the  structure  and  bonding  of  B4C  from 
previous  studies13,30,jl  is  as  follows.  Figure  la  shows  the  unit 
cell  of  B4C,  which  belongs  to  the  rhombohedral  space  group 
(R3m).  It  consists  of  one  icosahedron  per  cell  with  six  lattice 
vertices  (denoted  as  e  for  equatorial  site)  that  each  connect  to  a 
different  icosahedron  through  a  three  atom  chain  alon^  the 
(lll)r  direction  where  r  represents  rhombohedra.13’' 1  In 
addition,  each  icosahedron  makes  direct  connections  to  the 
other  six  neighboring  icosahedra  through  six  lattice  vertices 
denoted  as  p  for  polar  site.  The  chemical  formula  B4C  used  to 
refer  to  boron  carbide  is  understood  to  denote  a  range  of 
stoichiometries  because  the  locations  of  carbon  atoms  are  not 
determined  experimentally.  DFT  predicts  that  the  ground-state 
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Figure  1.  (a)  (BHCp)CBC  ground-state  structure  for  B4C  (rhombohedral  representation);  (b)  Si  microalloyed  structure,  (BuCp)— Si2;  and  (c) 
predicted  phonon  spectrum  for  (BnCp)— Si2.  Here  boron  is  green,  carbon  is  sienna,  and  silicon  is  blue. 


structure  is  (BuCp)CBC. 13,30  The  CBC  three-atom  chain  has 
two  two-center— two-electron  (2c— 2e)  sigma  bonds  from  the 
central  B  of  the  chain  to  the  chain  carbon  atoms,  while  each 
chain  carbon  also  makes  three  2c— 2e  bonds  to  boron  atoms  of 
three  adjacent  icosahedra.  Thus,  the  middle  boron  in  the  CBC 
chain  has  one  net  positive  formal  charge.  As  a  result,  the 
(BnCp)  cage  has  26  e  involved  in  intraicosahedral  bonding, 
stabilizing  the  icosahedral  structure  (Wade’s  rule).32'33 

To  examine  how  microalloying  might  modify  the  mechanical 
properties,  we  start  by  modifying  the  chain  structures 
connecting  icosahedra.  We  tried  many  systems,  but  the  most 
interesting  is  the  modified  structure  in  which  a  Si— Si  pair 
replaces  the  CBC  chain,  leading  to  the  chemical  formula 
(Bncp)  — Si2,  as  shown  in  Figure  lb.  The  ductility  improvement 
is  due  to  changes  in  the  deformation  mechanism  upon 
micro  alloying  with  Si. 

Several  alloyed  boron  materials  retaining  the  icosahedral 
motif  have  been  synthesized  experimentally,  including  B4C, 
B12P2,  and  B60.  These  are  formed  by  inserting  a  three-atom  or 
two-atom  chain  unit  into  six  of  inter  icosahedral  bonds  of  alpha 
boron.31  Our  studies  of  the  deformations  of  B4C  suggested  that 
the  interaction  between  chain  and  cage  plays  an  essential  role  in 
deconstruction  of  the  structures.30  Thus,  we  examined  how  the 
mechanical  properties  change  upon  modifying  the  CBC  chain 
in  (BnCp)CBC  structure.  In  particular,  we  replaced  the  C— B— 
C  chain  with  various  two-atom  pairs  or  single  atoms,  such  as 
Si— Si,  N— N,  S— S,  P— P,  CH2C,  and  S.  We  also  included  the 
high  energy  B4C  structure  (B12)CCC  (1.09  eV/unit  cell  higher 
in  energy)  and  the  boron-rich  boron  carbide  (B12)CBC. 

The  calculated  densities,  ductility  indexes,  and  hardness  are 
shown  in  Table  1.  We  find  that  (BnCp)— Si2  has  the  best 
ductility  index  of  1.59  compared  with  1.20  for  B4C  and  1.03  for 
a-boron.  It  also  has  a  low  density  of  2.454  g/cm3  (2.7%  lower 
than  the  2.523  g/cm3  of  (BnCp)CBC.  The  hardness  decreases 
slightly  to  27.8  GPa  compared  with  31.7  GPa  for  (BnCp)CBC. 
The  other  new  structures  with  improved  ductility  and  lower 
density  include  (B12)CCC,  (B12)CBC,  and  (BnCp)CH2C, 
which  have  B/G  equal  to  1.26,  1.37,  and  1.27,  respectively.  We 
observed  that  (BnCp)S  has  a  density  and  ductility,  similar  to 
(Bncp)  — Si2,  but  its  decreased  hardness  of  21.4  GPa  makes  it 
less  promising. 


Table  1.  Predicted  Density,  Ductility  Index,  and  Knoop 
Hardness  of  Materials  Modified  from  Boron  Carbide 


compound 

density 

(g/W) 

B/G  (ductility 
index) 

hardness 

(GPa) 

(BnCp)(C-B-C) 

2.523 

1.20 

31.7 

(BnCp)(Si— Si) 

2.454 

1.59 

27.8 

(B12)(C— B— C) 

2.439 

1.37 

30.2 

(B12)(C— C— C) 

2.470 

1.26 

41.4 

(BnCp)(NN) 

2.605 

1.11 

27.0 

(BuCp)S 

2.540 

1.53 

21.4 

(BuCp)(SS) 

2.651 

1.43 

23.6 

(BnCp)(P-P) 

2.642 

1.00 

24.0 

(BnCp)(C— H  H-C) 

2.404 

1.27 

33.9 

Because  the  B/G  criterion  indicates  that  (BnCp)— Si2  has  the 
best  ductility,  we  examined  the  stability  of  this  new  structure. 
The  PBE  predicted  structure  has  the  space  group  of  Cm  with 
primitive  cell  parameters  of  a  =  5.363  A,  c  =  5.099  A,  a  = 
109.207°,  and  y  =  70.622°.  The  predicted  phonon  spectrum 
(Figure  lc)  shows  that  the  Hessian  has  no  negative  eigenvalues, 
indicating  that  this  structure  is  dynamically  stable.  The  B/G 
criterion  is  only  an  empirical  rule  to  estimate  ductility  for 
metals.34  To  prove  that  the  new  structure  (BnCp)— Si2  has 
good  ductility  and  to  validate  the  B/G  criterion,  we  sheared  the 
system  along  the  amorphous  slip  system  of  (0111)/<1101> 
found  in  B4C.29,30  We  also  sheared  the  (B12)CCC  and 
(B12  )CBC  structures  expected  to  have  improved  ductility. 

Figure  2  shows  the  stress— strain  curves  for  these  three 
systems  along  with  (BnCp)CBC.  These  results  show  that  the 
failure  strain  is  0.348  for  (BnCp)CBC,  while  (B12)CCC  and 
(B12)CBC  lead  to  a  failure  strain  of  0.381  and  0.348.  The  ideal 
shear  strengths  for  (B12)CCC  and  (B12)CBC  are  51.0  and  45.8 
GPa  compared  with  39.0  GPa  for  (BnCp)CBC.  This  illustrates 
that  modifying  the  chain  structures  can  improve  the  ductility 
and  that  the  B/G  criterion  is  useful  for  predicting  changes  in 
ductility.  However,  current  methods  of  synthesis  cannot  control 
the  fractions  of  carbon  in  the  chains  and  icosahedra,  and 
current  materials  may  already  have  grains  or  regions  with  these 
more  favorable  materials.  Of  course,  the  grain  boundaries  may 
play  an  important  role  in  ductility,  a  factor  not  considered  here. 
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Shear  strain 

Figure  2.  Stress— strain  relation  of  various  structures  shearing  along 
the  (011l)/<1101>  amorphous  slip  system  of  B4C.  The  ground-state 
(B11Cp)CBC  fails  catastrophically  at  a  strain  of  0.348.  The  less  stable 
structure  (B12)CCC  and  the  boron-rich  structure  (B12)CBC  fail 
catastrophically  at  strains  of  0.348  and  0.381.  The  ductile  structure 
(BnCp)  — Si2  shows  structural  changes  at  strains  of  0.397  and  0.741 
without  brittle  failure.  The  filled  circles  represent  the  key  structural 
steps  examined  in  Figures  3  and  4 


We  next  consider  that  a  system  such  as  (BnCp)— Si2.  The 
predicted  stress— strain  curve  for  (BnCp)— Si2  (Figure  2)  shows 
a  plastic  ductile  behavior.  That  is,  it  is  elastic  until  a  shear  of 
0.245,  but  instead  of  failing  by  amorphous  shear  as  in  B4C,  we 
see  the  stress  continuous  increase  after  strains  of  0.397  and 
then  decrease  at  the  strain  of  0.715,  about  twice  that  for  B4C. 
The  icosahedra  do  not  break  during  the  shear  to  the  large  strain 
of  0.802. 

To  understand  these  structural  transformation  mechanisms, 
we  examined  its  relation  to  the  bonding  changes  for  (BnCp)— 
Si2;  Figures  3  and  4  show  the  structures  and  the  isoelectron 
surface  (at  0.9)  from  the  electron  localization  function  (ELF) 
analysis  at  various  critical  strains.  The  ELF,  which  ranges  from  0 
to  1,  enables  an  effective  and  reliable  analysis  of  covalent 
bonding  and  lone  pair  formation.35,36 

Figure  3  shows  the  first  structural  transition  to  another  stable 
conformation.  Figure  3a  shows  the  intact  structure  where  the 


icosahedron  has  25  electrons  for  intraicosahedral  bonding,  one 
less  than  Wade’s  rule  because  each  Si  makes  four  2c— 2e  bonds. 
As  the  strain  increases  to  0.245  that  corresponds  to  maximum 
stress  of  30.4  GPa,  the  B— C  bonds  between  icosahedra  increase 
by  20%  to  1.997  A,  as  shown  in  Figure  3b.  Continuous  shear  to 
0.263  strain  leads  to  the  stress  drops  through  breaking  B— C 
bond  between  cages,  as  shown  in  Figure  3c.  The  electron 
counting  scheme  at  0.263  strain  changes  and  the  cage  has  27 
electrons.  As  the  strain  increases  to  0.314  (Figure  3d),  the  Si- 
Si  bond  in  the  chain  breaks  and  the  Si  atoms  rebond  to  the 
unbonded  B  and  C  atoms  in  the  cage.  Now  the  cage  has  25 
electrons  again.  As  the  strain  increase  to  0.397  (Figure  3e), 
another  stable  structure  forms  where  one  Si  connects  to  three 
boron  atoms  in  cages  and  the  other  Si  connects  to  two  boron 
and  one  carbon  atoms  in  cages.  This  new  stable  structure  is 
slightly  lower  in  energy  by  7.9  kcal/mol  per  icosahedron. 

Figure  4  shows  the  second  structural  transition  where  the 
new  structure  recovers  to  its  original  structure.  Figure  4a  shows 
the  structure  that  corresponds  to  the  second  maximum  stress  at 
0.689  strain.  The  B— B  bonds  between  cages  are  stretched  to 
2.18  A.  As  the  shear  increases  to  0.715,  the  Si— Si  chain  bonds 
and  B— B  bonds  between  cages  break  and  the  Si  rebonds  to  the 
unbonded  boron  atom  in  the  cage,  as  shown  in  Figure  4b. 
Continuous  shear  to  0.741  leads  to  the  reformation  of  Si— Si 
chain  bond  and  B— B  bond  between  cages,  as  shown  in  Figure 
4c, d.  It  is  worth  noticing  that  the  electron  counting  scheme 
does  not  changes  during  the  second  structural  transition. 

There  is,  of  course,  a  practical  problem  in  synthesizing  this 
(Bn  Cp)— Si2  material.  How  do  we  ensure  that  the  Si  does  not 
go  inside  the  icosahedron  and  the  C  and  B  do  not  go  into  the 
chain?  Experimental  synthesis  from  various  starting  materials 
involves  complex  transformations  well  beyond  the  QM 
methods  used  here.  However,  the  theory  can  help  by  predicting 
signatures  for  the  atoms  being  at  various  locations.  For 
example,  the  QM  calculated  NMR  chemical  shifts  for 
(BnCp)-Si2  are  86.4  ppm  for  Cp  and  457.3  and  475.6  ppm 
for  Si,  whereas  for  (BnSi)-CSi  they  are  36.2  ppm  for  Cchain, 
464.2  ppm  for  Sichain,  and  814.8  ppm  for  Sicage,  respectively. 
Thus,  by  measuring  XPS,  XANES,  NMR,  Raman,  and  IR  for 


(d)  Strain  =  0.314  (e)  Strain  =  0.397 

Figure  3.  (BuCp)  — Si2  structures  and  ELF  of  the  first  structural  transition  shearing  along  (011l)/<1101>  amorphous  slip  system  of  B4C:  (a)  the 
intact  structure;  (b)  the  structure  corresponds  to  the  maximum  stress  where  the  B— C  bond  distance  between  icosahedra  increases  by  20%  to  1.997 
A,  point  A  in  Figure  2;  (c)  the  structure  where  the  B— C  bonds  between  icosahedra  break,  point  B  in  Figure  2;  (d)  the  structure  where  the  Si— Si 
bonds  in  the  chain  break  and  the  Si  atoms  rebond  to  the  unbonded  B,  C  atoms  from  previous  breaking  B— C  bond,  point  C  in  Figure  2;  and  (e)  the 
transformed  structure  with  new  formed  Si— Si  bond  and  B— B  bond  between  icosahedra,  point  D  in  Figure  2 
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(c)  Strain  =  0.741  (d)  Strain  =  0.741 


Figure  4.  (BnCp)  — Si2  structures  and  ELF  of  the  second  structural  transition  shearing  along  (011l)/<1101>  amorphous  slip  system  of  B4C:  (a)  the 
structure  corresponds  to  the  second  maximum  stress  and  the  B— B  bond  between  cages  (red  dashed  line)  is  stretched  to  2.18  A,  point  E  in  Figure  2; 
(b)  the  structure  where  the  Si— Si  chain  bonds  and  B— B  cage  bonds  break  and  Si  rebonds  to  the  unbonded  boron  atoms  in  the  cage,  point  F  in 
Figure  2;  (c)  the  second  structural  transition  where  the  Si— Si  chain  bond  and  B— B  cage  bonds  reform,  point  G  in  Figure  2;  and  (d)  the  same 
structure  as  (c),  but  viewed  along  “a”  direction  that  clearly  shows  the  reformation  of  Si— Si  chain  bond. 


samples  prepared  under  various  synthesis  conditions,  we  might 
be  able  to  select  conditions  that  favor  (BnCp)— Si2  by 
comparing  to  the  predicted  signatures  for  each  distribution  in 
the  experiment. 

In  summary,  we  show  that  mechanical  properties  can  change 
dramatically  as  we  modify  the  CBC  chain  in  the  (BuCp)CBC 
structure.  In  particular,  replacing  the  CBC  chain  with  the  Si— Si 
pair  dramatically  increases  the  ductility  of  the  boron  carbide. 
Moreover,  the  density  decreases  by  2.7%,  while  the  hardness 
decreases  by  only  12%. 

■  COMPUTATIONAL  METHODS 

All  DFT  calculations  were  performed  using  the  VASP 
package,37  39  with  Perdew— Burke— Ernzerhof  (PBE)  func¬ 
tional  and  the  projector-augmented  wave  method  to  account 
for  the  core— valence  interactions.40  We  found  that  a  kinetic 
energy  cutoff  of  600  eV  for  the  plane-wave  expansions  gives 
excellent  convergence  of  the  total  energies,  energy  differences, 
and  structural  parameters.  Reciprocal  space  was  sampled  using 
the  r -centered  Monkhorst— Pack  scheme  with  a  resolution  less 
than  In  X  1/40  A-1.  The  convergence  criteria  were  set  to  1  X 
10-s  eV  energy  difference  for  solving  for  the  electronic  wave 
function  and  1  X  10-3  eV/A  force  for  geometry  optimization. 

To  establish  the  stability  of  the  predicted  crystal  structure,  we 
calculated  the  Hessian  and  phonon  modes  using  the  supercell 
approach  with  finite  displacements,  as  implemented  in 
PHONPY  code.41 

To  make  a  preliminary  estimate  of  ductility,  we  used  the 
Paugh  criterion'4  that  high  values  of  B/G  (particularly  B/G  > 
1.75)  are  more  ductile,  where  B  is  the  bulk  modulus  and  G  is 
the  shear  modulus.  We  first  computed  the  elastic  constant  C ^ 
from  the  stress— strain  relation  as  a  function  of  various  cell 
distortions  from  the  equilibrium  lattice  configuration.42  Then, 
the  inverse  of  the  elasticity  tensor  was  derived  from  Sl;  = 
( Cy)  1 .  Finally,  we  used  the  Voigt— Reuss— Hill  (VRH) 
approximation43  to  calculate  the  isotropic  polycrystalline  elastic 
moduli  from  the  corresponding  single-crystal  elastic  constants. 
The  VRH  average  is  the  arithmetic  mean  of  the  Voigt  and 


Reuss  bounds.  For  Voigt  average,  the  bulk  modulus,  Bv,  and 
shear  modulus,  Gv,  are  expressed  as 

Bv  =  (Cn  +  C22  +  C33  +  2C12  +  2C13  +  2C23)/9 

and 

Gv  =  (Cn  +  C22  +  C33  —  C12  —  C13  —  C23  +  3C44 
+  3C55  +  3CJ/15 

For  Reuss  average,  the  bulk  modulus,  Br,  and  shear  modulus, 
Gr,  are  expressed  as 

=  l/(Sn  +  S22  +  S33  +  2S12  +  2S13  +  2  S23) 

and 

Gr  =  15/ (4S, ,  +  4S22  4-  4S33  —  4S12  —  4S13  —  4S23 
+  3S44  +  3S55  +  3S66) 

We  also  computed  the  Knoop  hardness  based  on  the 
electronegativity,  covalent  radii  of  the  constituent  atoms,  and 
the  bond  lengths  in  the  structure.44,45 

To  examine  the  shear  failure  of  the  modified  materials,  we 
imposed  the  strain  for  a  particular  shear  plane,  while  allowing 
full  structure  relaxation  of  the  atoms  and  of  the  other  five  strain 
components.46  The  residual  stresses  after  relaxing  these  other 
strain  directions  were  <0.5  GPa.  For  these  shear  simulations,  we 
used  supercells  varying  from  112  atoms  to  120  atoms,  and  the 
Monkhorst— Pack  grid  (3  X  3  X  3)  in  the  k-space  was  used. 
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